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Abstract 

The Imaginary Time Step (ITS) method is apphed to solve the Dirac equation with the nonlocal 
potential in coordinate space by the ITS evolution for the corresponding Schrodinger-like equation 
for the upper component. It is demonstrated that the ITS evolution can be equivalently performed 
for the Schrodinger-like equation with or without localization. The latter algorithm is recommended 
in the application for the reason of simplicity and efficiency. The feasibility and reliability of this 
algorithm are also illustrated by taking the nucleus ^^O as an example, where the same results as 
the shooting method for the Dirac equation with localized effective potentials are obtained. 
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With the operating of the worldwide new Radioactive Ion Beam (RIB) facilities and 
the developments in the detection techniques, a lot of novel aspects of nuclear structure and 
entirely unexpected features have been found in the exploration of "exotic nuclei" [l, 2, ^, 4|. 
The extreme neutron richness of exotic nuclei and the physics related to the low density in 
the tails of nuclear matter distributions provide both demand and challenge to solve complex 
many-body problem in coordinate space. 

As one of the best candidates for the description of exotic nuclei, the Relativistic Mean 
Field (RMF) approach Slhas achieved lots of success in describing many nuclear phenomena 
during the past years 6, 7, ^. In particular, the Relativistic Continuum Hartree-Bogoliubov 
(RCHB) theory 9|] provides a fully self-consistent treatment of pairing correlations in the 
presence of the continuum and thus gives a reliable description of nuclei far away from the line 
of /3-stability, which suggests a new mechanism for the formation of the halo phenomena 10|. 



Within the RMF approach, however, the Fock terms of the energy density functional are 
neglected under the Hartree approximation. The recent development of density-dependent 
relativistic Hartree- Fock (DDRHF) theory ll| resumes the Fock terms and obtains quantita- 



tively comparable precision with the RMF theory for nuclear matter and finite nuclei H 



12|. 



Moreover, within the DDRHF theory, the importance of the Fock terms is evidenced by the 
improvement on the descriptions of the nuclear shell structures 12| and their evolutions 13| , 
the excellent reproduction of the spin-isospin resonance based on the self-consistent random 
phase approximation [l^, as well as the influences on isospin properties of nuclear matter 



and neutron stars at high densities [15|]. More recently, as an extension of the DDRHF the- 
ory, the relativistic Hartree-Fock-Bogoliubov theory with density dependent meson-nucleon 
couplings (DDRHFB) is developed for the description of the exotic nuclei jl^. 

For the future study, one of the natural extensions of the RCHB or DDRHFB theory is 
the exploration of deformed exotic nuclei, which will provide the decisive conclusion on the 
existence of deformed halo. Efforts along this line have been made in the past several years. 
Due to the difficulty in solving the coupled channel differential equations for deformed system 
in coordinate space the method for expansion in Woods-Saxon basis was proposed 18|, 
which becomes very time consuming for the heavy system. 

The Imaginary Time Step (ITS) method 19| is another effective approach to solve the 
nonrelativistic problems in coordinate space and has achieved lots of success in the conven- 
tional mean field approach [2^ . More importantly, its extension from spherical to deformed 
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systems is straightforward. It is therefore worthwhile to adopt this method for the relativis- 
tic system. Naively, one may assume that the ITS method will result in a great disaster due 
to the Dirac sea. Fortunately, it has been demonstrated that this disaster can be avoided 
by the ITS evolution for the Schrodinger-like or charge-conjugate Schrodinger-like equation 



for the solution respectively in the Fermi or Dirac sea [21 1. 

For the application of the ITS method in the relativistic Hartree-Fock theory, one has to 
solve the Dirac equation with nonlocal potential due to the Fock terms. The corresponding 
Dirac equation will be a set of coupled integro-differential equations, which cannot be solved 
by the conventional shooting method directly. 

In this letter, the ITS method will be applied to solve the Dirac equation with the nonlocal 
potential in coordinate space following the success in solving the Dirac equation with the 
local potential by the ITS evolution for the corresponding Schrodinger-like equation for the 
upper component 2l|]. It will be demonstrated that the ITS evolution can be equivalently 
performed for the Schrodinger-like equation with or without localization, and the latter 
algorithm is recommended for the reason of simplicity and efficiency. 

In the relativistic Hartree-Fock approach, the main issue is to solve the Dirac equation 

{a • p + /3 [M + Sir)] + V{r)} v.„(r) + J dr'U{r, r')ifa{r') = s,ifa{r), (1) 

with the local (scalar and vector) potentials S{r) and V{r), and the nonlocal potential 
U{r,r'). For clarity and simplicity, V(r) ± S{r) are assumed to be the spherical Woods- 



Saxon potentials similar as in Ref. 18|], and thus the Dirac spinor takes the form. 



Mr) = t ( ^"^l I X. (g.), la + l = 2j., (2) 



-Ga{r)y]:.^jn) 

with the isospinor Xi(Q'a) aiid the spherical spinor y'-^.[Q). The single-particle state is 
labeled by the quantum numbers {ria, /q, k^, ja} of the upper component, where = ±|ja + 
1/2| for =j,±l/2^ 



Similar as in Ref. 22|, the nonlocal potential U{r,r') is expanded as. 



ui{r,r') 2/ + 1 
rr' An 

where Pi is the Legendre polynomial and u is the angle between r and r'. More explicitly, 



f/r,r => Piicosu), 3) 

rr An 

I 
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if r') takes the same form as in Ref. 22l |. 



f/(r,rO = ^53$0^_(r.-nW^ (4) 



the expansion coefficient ui{r,r') can be obtained as 

2\/rr' I -T—kiH. ( Irr' 



where J;^i is the modified Bessel function of the first kind. In practice, the central part u{r) 
could take the shape of Woods-Saxon type, 

with the parameters, Vq, ro, a and the nonlocality parameter 7. 

With the expansion as in Equation ([3]), the Dirac equation with nonlocal potential can 
be reduced as 

V + S + M -i + f \ f F^{r)\fx^{r)\f F^ir)\ 
+ £ + V-{S + M) ) \Ga{r) ) \Ya{r) ) 1^ G«(r) j ' 

where the nonlocal terms are 

Xa{r) = J dr'Ux{r,r')F,{r'), Ux{r,r') = Ui^{r,r'), (8a) 

Ya{r) = J dr'Uy{r,r')Ga{r'), Uyir^r') = ui,{r,r'). (8b) 

The coupled integro-differential equations in Equation([7]) cannot be directly solved by the 
conventional shooting (Runge-Kutta) method unless they are transformed to homogeneous 
differential equations by localization 



231 . For example, one can rewrite Equation ([7]) as 



dr 

l_lV„ T/_ 

dr r 

with the localized effective potentials. 





















\ Ga{r) j 




\ Ga{r) j 



(9) 



Y - Y - (lOh) 



aV/'aV/ a' 



then Equation([9]) could be solved iteratively from an initial guess of the wave function by 
the shooting method. 



For the ITS method, as demonstrated in Ref. 2l(], the evolution can be performed by using 
the effective single-particle Hamiltonian for the upper component h^g . With the relation 
between the upper and lower components in Equation 

^- = Tr + + Ya,F.F^ , where Ml = M - {V - S) - Ya,G. + (11) 

Ml \ dr r J 

the Schrodinger-like equation for the upper components is, Kb Fa = SaFa, with 



^off. Fa 



1 d^Fa 



+ 



+ 
+ 



Ml dr^ 

1 Kajhia + l 

Mr r2 



^ (Y X 1+ ^ 



dFa 



dr 



+ 



1 Kg / 1 dML 
Ml r \Ml dr 



Xa. 



Ga 



1 dYa,F^ 1 dML 

Ml dr Ml dr '"'^^ 



1 



Mr 



Xa^Gja^F. + (\/ + 5 + M + 



(12) 



On the other hand, one could also obtain the relation between the upper and lower 
components from Equation ([7]), 

1 



Gn 



M^ 



dFa , Ka „ 
\ 

dr r 



Fj, where M+ = M - (^ - 5) + 



(13) 



and the corresponding Schrodinger-like equation for the upper component gives 



Ks.Fa 



1 d^Fa 



+ 



1 dM<dFn 



M+ dr^ Ml dr dr 



+ 



/T. ^ , ,N 1 dM+ Ka 1 KaiHa + l] 

{V + S + M) + , — + ^ 



Ml dr r 



M~ 



1 dYa 



+ 



1 dM^ 
Ml dr 



'-Ya + 



1 K„ 



-Ya + Xa. 



Fa 
(14) 



M+ dr ' Ml dr ' M+ r 

Although the right-hand sides are different. Equation ( fT4|) is exactly the same as Equa- 
tion (fT2l) . Moreover, one should note that the ITS evolution involves only the operation 
h^gFa- Therefore, the nonlocalized ITS evolution for Equation (|T4|) is equivalent to the 
localized one for Equation f|T2l) . For the reason of simplicity and efficiency, the nonlocalized 
ITS evolution for Equation (fT^ will be used. 

In the numerical calculations, the ITS evolution starts from the initial orthogonal wave 
functions {^p^a^} with spherical Bessel function for the upper components and zero for the 
lower components. Together with the given nonlocal potentials (l5l), the initial X^\r) and 
Yj^\r) are constructed according to Eqs. (IHaj) and ( I8bl) . Then the operation h^s.Fa in 
Equation (fT4l) is obtained and applied to evolve the single-particle wave functions, which is 
carried out in coordinate space within a spherical box [0, R] with the box size R = 20 fm, 



mesh size dr = 0.1 fm and a typical time step At = 10~^^ s. The set of wave functions 
thus obtained is then orthogonahzed by the standard Gram-Schmidt procedure to provide 
a new set of orthogonal single-particle wave functions which will be used for the next ITS 
evolution. This process is repeated till the final convergence. 

Take the state v\s\i2 in ^^O as an example, the nonlocal potentials Ux{r, r') and f/y(r, r') 
defined in Eqs. (l8a|) and (18bp are shown in Fig. [1] The detailed parameters for the nonlocal 
potentials chosen tentatively to give a reasonable single-particle spectrum in ^^O can be 
found in the corresponding caption. The nonlocal potentials Ux{f-,r') and Uyir^r') are 
state-dependent with different orbital angular momentum quantum numbers la and l'^ for 
the upper and lower components, respectively. 

The ITS evolutions of the single-particle wave function and the nonlocal terms X{r), Y{r) 
for vlsi/2 in ^^O are illustrated in Fig. [2l The initial single-particle wave function which 
is assumed to be spherical Bessel function for the upper component and zero for the lower 
component, together with the corresponding X{r) and Y{r) are denoted by the thin solid 
line, which are followed by the results at 2 x 10'^ (dotted lines) and 4 x 10^ (dash-dotted lines) 
iterations. Finally, the evolution results converge to the same ones obtained by the shooting 
method. The results from shooting method are obtained by solving the Dirac equation with 
the localized effective potentials from the converged X{r) and Y{r) in the ITS evolution. 

Similarly, the neutron single-particle spectrum of ^^O can be obtained by the ITS method 
as shown in Fig. [3l where the the evolutions of the corresponding single-particle energies are 
also illustrated. It can be seen that different levels converge at different speeds. The deeper 
the level lies, the faster the evolution converges. Most levels converge monotonously except 
the level 2s 1/2, which evolves with fluctuation at the beginning due to the orthonormalization 
with respect to the level lsi/2. 

In summary, the ITS method is applied for solving the Dirac equation with nonlocal po- 
tential. Similar as in Ref. 2l|], the evolution is performed for the corresponding Schrodinger- 
like equation for the upper component. It is demonstrated that the ITS evolution can be 
equivalently performed for the Schrodinger-like equation with or without localization. The 
nonlocalized ITS evolution is much more simple and efficient, and thus is recommended for 
future application. These conclusions are further supported by the numerical solution of the 
Dirac equation with nonlocal potential reasonably chosen for ^^O. The present algorithm 
provides the same results as the shooting method for the Dirac equation with the localized 
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effective potentials. The investigation here demonstrates the possibihty of the ITS method 
in the relativistic Hartree-Fock theory. 
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FIG. 2: (color online) The evolutions of the single-particle wave functions F{r), G{r) and the 
nonlocal terms X{r)^ for ^^^^1/2 iii nucleus ^^O. The evolutions are illustrated at the initial 
(thin solid lines), 2 x 10^ (dotted lines), 4 x 10^ (dash-dotted lines) iterations. The converged 
results (thick solid lines) of the ITS evolution are compared with the results by shooting method 
(dashed lines). 
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FIG. 3: (color online) The evolutions of the single-particle energies and the final convergent 
neutron single-particle spectrum of nucleus ^^O obtained by the ITS method. 
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